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1 Introduction 

The Driver integration by parts formula [12] and the Bismut derivative formula [4] are two 
quite useful tools in various aspects of stochastic analysis. Let V be the gradient operator 
and Pt stand for the diffusion semigroup. The above two formulas allow us to estimate the 
commutator VP( — P*V, which plays a key role in the study of flow properties [16j . On the other 
hand, [39] showed that, in general the integration by parts formula is more complicated and 
harder to obtain than the derivative formula. Based on martingale method, coupling argument 
or Malliavin calculus, the derivative formula has been widely studied and applied in various 
fields, such as heat kernel estimates, strong Feller property and functional inequalities, see 
[HI [37J EH H3] and references therein. Whereas, in [39], based upon a new coupling argument, 
the integration by parts formulae are derived and applied to various models including degenerate 
diffusion process, delayed SDEs and semi-linear SPDEs. 

Recently, transportation cost inequality has been widely studied. Let (E, d) be a metric 
space equipped with cr-algebra 33 such that d(-,-) is 38 x 38 measurable. For any p > 1 and 
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two probability measures // and v on (E, 3$), the L p -Wasserstein distance induced by d between 
these two probability measures is defined by 

W$M= inf (If d{x,yf-K{dx,dy) 
ttGV(m,u) \JeJe 

where ^(/i, v) denotes the set of all coupling of fj, and v. In 1996, Talagrand [HI] proved the 
following transportation cost inequality for the standard Gaussian measure [i on M. d : 

WiUix^f < 2/i(/ log/), / > 0,/i(/) = 1, 

where d(x,y) = \x — y\. In general, we call that the probability measure ji satisfies the L p - 
transportation cost inequality on (E,d), if there exists a constant C(> 0) such that for any 
probability measure v, 

W*M< (1.1) 
where H(z/|/i), the relative entropy of v with respect to /t, is given by 

[ +oo, else. 

For simplicity, we write [i G T p (C|d) for (jl.ip . In the past decades, the work of Talagrand 
has been generalized to various different stochastic processes, see, for instance, \27\ [5] for the 
Hamilton- Jacobi equation, [171 EB] on abstract Wiener space, [15] on loop groups, [TT 1 [34" 1 |35 | [36] 
for diffusion processes, [33] for SDEs of pure jumps, [21] for SDEs driven by both Brownian 
motion and jump process, [3] for neutral functional SDEs, [30] for SDEs driven by fractional 
Brownian motion. 

In this article, we are interested in the stochastic differential equations driven by a fractional 
Brownian motion. It is well known that the main difficult point raised by the fractional Brownian 
motion is that it is not Markovian process and semimartingale, so the ltd approach to setup 
a stochastic integral with respect to the fractional Brownian motion is not valid. Now there 
exist numerous attempts to define a stochastic integral with respect to the fractional Brownian 
motion and moreover, many works to discuss the stochastic differential equations driven by a 
fractional Brownian motion. We briefly present some results. Based on a fractional integration 
by parts formula |42j . Nualart and Ra§canu |25] established the existence and uniqueness result 
with H > \. By using the theory of rough path analysis introduced in [20] , Coutin and Qian [7] 
proved an existence and uniqueness result with Hurst parameter H € ( 5)- -F° r the regularity 
results about the law of the solution, one can see [191 123"! [26] and references therein. 

The equation driven by a fractional Brownian motion we are to deal with is of Volterra type, 
which is originally discussed by Coutin and Decreusefond [6]. In [6], they studied existence, 
uniqueness and regularity of solution. In this paper, by using coupling argument and the Gir- 
sanov transform for fractional Brownian motion, we first prove the Driver integration by parts 
formula for aimed equation. As an important application, we give an alternative proof for [61 
Corollary 4.1]. Secondly, the Bismut derivative formula is established, using the techniques of 
Malliavin calculus. As applications, gradient estimate, the Harnack inequality and the strong 
Feller property are derived. Finally, we obtain the Talagrand type transportation cost inequal- 
ities for the law of the solution of aimed equation on the path space with respect to both the 
uniform metric and L 2 -metric. 




2 



The paper is organized as follows. In section 2, we give some preliminaries on fractional 
Brownian motion. In section 3, we investigate the Driver integration by parts formula, while 
in section 4, the Bismut derivative formula is discussed. Finally, section 5 is devoted to the 
transportation cost inequalities. 

2 Preliminaries 

Let B H = {B^ 1 ,t E [0, T]} be a fractional Brownian motion with Hurst parameter H S (0, 1) 
defined on a complete probability space (CI, & , P). Namely, B H is a centered Gauss process with 
the covariance function 

E(B t H Bf ) = R H (t, s) := \ (t 2H + s 2H - \t - s\ 2H ) . 

When H = ±, the process B^ is the usual Brownian motion. By the above covariance function 
and the Kolmogorov criterion, we know that B H have (H — e)-order Holder continuous paths 
for all e > 0. Furthermore, B H has stationary increments and is self-similar with Hurst index 
H. 

From [10], it is known that the covariance kernel Ru(i, s) admits the following representation: 

ptAs 

R H (t,s)= K H (t,r)K H (s,r)dr, 
Jo 

where Kh(-, •) is a square integrable kernel given by 

K H (t, s) = t(h + ±\ (t- s) H ~^F (h - I \ - H, H + i 1 - ^ , 

in which F(-, ■, •, •) is the Gauss hyper geometric function (for details, see |22j). 

Again by [10], the operator K H : L 2 ([0, T];R) — > I^ 1/2 (L 2 ([0, T];R)) associated with the kernel 

Kh(-, •) is defined as follows 

(K H f)(t) := / K H (t,s)f(s)ds, 
Jo 

where I 0+ is the (H + l/2)-order left fractional Riemann-Liouville integral operator on [0, T]. 
It is an isomorphism and for each / G L 2 ([0,T];M), 

(K H f)(s) = I 2 *sV 2 - H iy;- H s H -lj, H < 1/2, (2.1) 
(K H f)(s) = ll + s H -V 2 I»- ll2 sV 2 - H f, H > 1/2. (2.2) 

Hence, for any h £ I^ 1 ^ 2 (L 2 ([0, T]; R)), the inverse operator K^ 1 can be written as 

(K^h)(a) = s^D^s^h', H > 1/2, 
(K^h)(s) = s^ H DH 2 - H s H -y 2 D 2 »h, H < 1/2, 

where Dq + is the a-order left-sided Riemann-Liouville derivative operator, a € (0,1). 
In particular, when h is absolutely continuous, it holds 

(Kjj l h)(s) = s H ^ 2 I 1 / + 2 - H s 1 ^ H h', H < 1/2. 
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For more details about the deterministic fractional calculus, one can refer to |29j. 

We assume that Q is the canonical probability space Cb([0, T]; R), the set of continuous 
functions, null at time 0, equipped with the Borel cr-algebra and P is the law of the fractional 
Brownian motion. The canonical filtration is J^t = <^{B^ : < s < i] V Af, where Af is the set 
of the P-null sets. According to |10| Theorem 3.3], the Cameron-Martin space of the fractional 

Brownian motion, denoted by %, is equal to / ( ^ +1/ ' 2 (L 2 ([0, T]; R)), i.e., for any h G %, it can be 
represented as h(t) = Kffh(t), where the function h belongs to L 2 ([0, T]; R). The scalar product 
on H is defined by 

(h,g)n ■= (KH 1 h,KH 1 9)L a ([o,T\;TSL), VKg e U. 

As a consequence, (f2,%,P) is an abstract Wiener space in the sense of Gross. Furthermore, let 
Q* denote the strong topological dual of f2, then there hold 

n* l 2 ([o,t];R) ^n^n 

and 

R H = K H oK* H , 

where we identify the operator Rjj and its kernel. 

Next we summarize some basic results of Malliavin calculus associated with the fractional 
Brownian motion, and we refer to [ID], [24] and |32] for a comprehensive presentation. 

Let S denote the set of smooth and cylindrical random variables of the form: 

F = f((h,0j),--;(l n ,0j)) 

where n > 1, / € C^°(W l ), the set of / and all its partial derivatives are bounded, U € 0*, 1 < 
i < n. The Malliavin derivative of F, denoted by DhF, is defined as the "H-valued random 
variable 

D h F(cj) = V %f-({h,u), ■ ■ (l n ,io))R H (k). 

For any k € N, denote by the iteration of Djj. For any p > 1 and k G N, we define the 
Sobolev space U>^f as the completion of S with respect to the norm: 

k 

\\F\\l >p :=E\F\r + EY,\\DhF\F Hm . 
i=i 

The divergence operator Sh, also called the Skorohod integral, is defined by using the duality 
relationship. More precisely, the domain Dom p <5f/ is the set of process u such that 

\E{D H F,u) H \ < C(E|F| 9 )i 

for all F G where q satisfies 1/p + 1/q = 1 and C is some constant depending on n. 
If u G Dom p (5fl-, then 5hu is defined by 

E(F^n) =E{D H F,u) n . 
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It is well known that, in the case of the Brownian motion (H = 1/2), the Skorohod integral 
is an extension of the ltd integral. So, this motivates us to use the divergence operator to define 
a stochastic integral with respect to the fractional Brownian motion. That is, 

f u t S H B^ :=5 H (K H u), 
Jo 

where the process Kjju G Dom<5# := U p >iDom p J^ (see e.g. [TO] and |6j). According to [TOl 
Theorem 4.8], we have the following Levy-Hida representation: 

W:=( I I [0A (s)8 h bA = (5 H (K H l m )) <t<T 

\Jo / 0<t<T 

is a standard Brownian motion whose filtration is equal to {^,0 < t < T} and moreover, for 
any u G L^([0, T] x fi), the set of square integrable and adapted processes, it holds 

u s 8 H B?= [ u s dW s ,\/t G [0, T]. 
Jo 

In particular, for each t G [0, T], taking u = Kjj(t, •), then we get B^ = J* Kn{t, s)dW s . 

In the paper, we are concerned with a R- valued equation driven by a fractional Brownian 
motion of the form: 

X t = x+ [ K H (t,s)b(s,X s )ds+ [ K H (t,s)a{s,X s )dW s . (2.3) 
Jo Jo 

Note that, when a = C, then the third term of the right-hand side of (|2.3p is equal to CB^ 1 . 
Hence, the factor KH(t,s) in the noise term is necessary to make the equation sense. While 
Kjj(t,s) in the drift term is only to symmetrize b and a. Set Hq = \H — l/2\, Ajj = {p > 1 : 
pH < 1} and for every p G A H , put k v = (1 - pH )-\ L K + P = U q>Kp L*([0, T];R). 

Definition 2.1 A M-valude process (^t)te[o,Tl IS called a solution of (|2.3p . if it is adapted such 
that E\X t \ 2 G L+ 2 , and §Z3§ is satisfied dP x dt a.s. 

Remark 2.2 From {21 Theorem 3.1 and Theorem 3.2], we know that, ifb and a are Lipschitz 
continuous for the second variable uniformly with respect to their first variable, and there exist 
x ,y G R such that b(-,x ) G L+ 1 ([0, T]; R), <r(-, y ) G L^ 2 ([0,T];R), then ([23]) has a unique 
solution. Furthermore, for allp G Ah, sup tg [ T ] E|X t | p < oo. If cr is bounded, then X has almost 
surely continuous trajectories. 

Define P t f(x) := E/(Xf), t G [0,T], / G 3S h {R), where is the solution to ([23]) with 
Xq = x and ^(R) denotes the set of all bounded measurable functions on R. Besides, we 
denote by C^(R) the set of all bounded continuous differentiable functions. In the remainder of 
the paper, we will establish the integration by parts formula and the Bismut derivative formula 
for Pt, and moreover obtain the Talagrand type transportation cost inequalities for the law of 
the solution of (|2.3|) on the path space. 
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3 Integration by parts formula 

This section is devoted to the equation (|2.3p with additive noise, i.e. 

X t = x+ f K H {t,s)b(s,X s )ds+ f K H {t,s)a(s)dW s . (3.1) 
jo Jo 

We aim to establish integration by parts formula by the method of coupling and Girsanov 
transformation. As an application, we give an alternative proof for [61 Corollary 4.1], in which 
the absolute continuity of the law of the solution is discussed. 

To start with, let us give some conditions of the coefficients b and a: (HI) 

(i) b is continuously differentiable w.r.t. the second variable and there exist positive constants 
K\ and K 2 such that 

\db(t,-)(x)\ < Ki,\a(t)- l \ < K 2 , V< G [0,T],x G R; 

(ii) there exist x G R such that 6(-,x ) G L+ 1 ([0, T]; R), <r(-) G L+ 2 ([0, T]; R). 
Theorem 3.1 Let T > and y G R be fixed. Assume that (HI) holds. 

(1) For each f G C£ (R), there holds the integration by parts formula 



P T (Vyf)(x)=E 



f(X£) J ais)' 1 (c H s l 2~ H - sdb(s, .)(*,)) |dW s 



where C# is a positive constant given in the proof below 
As a consequence, for each a > and positive f G C\ 



|iMV„/)| < a [P T (/ log f)-(Prf) (log P T /)] 

a V {2-2H)T 2H h-2H 3 



(2) For each non-negative / G ^(R), there holds the shift Harnack inequality 

pKl ( C% + 4C„K l T ,_„ K\T ■ y2 



(Pr/F < <«■{/(» 4- •)« ex P [p _ j V(2 _ 2H)T2H 5 _ 2ff - - 3 
(3) For each positive f G =^(R), there holds the shift log-Harnack inequality 
Pt log / < log P T {f(y + •)} + K\ ( - ^-^ + ^M"* + 



(2-2H)T 2H 5-2H 

Proof. Obviously, by (HI), it follows from Remark 12.21 that (13. lj) has a unique solution. On 
the other hand, for any e G [0,1], let X\ solve the equation 

X t £ = x+ f K H (t,s)b(s,X s )ds + f K H (t,s)a(s)dW s + %y, t e [0,T]. (3.2) 
Jo Jo 1 
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It is easy to see that X\ = X t + ^y, t G [0, T\. In particular, X^, = Xt + ey- 
Next, let 

( r(2ff)r(l/2-H) (l_rr\ rr . 1. 
B(2-2H,2H) \2 12 J ' n ^ 2 ' 

1, # 



rrg-i/2) 



l. 

2' 

H>i. 



K. B{2-2H,H~l/2) ' 

It follows from {HQ and (2.2) that, for each H G (0, 1), K H (C H x 1/2 ~ H )(t) = t. Therefore, we 
can reformulate (|3.2|) as 

Xl = x+ I K H {t,s)b(s,X e s )ds + [ K H {t,s)a(s)dW e s , te[0,T\, 
Jo Jo 

where W e s =W S + j s Q cr(r)- 1 (b(r, X r ) - b(r, X^ + ^pM"^ dr, s G [0, T]. 
Let 



&(r) = 6(r,X r ) - b(r,X*) + 



and 



R e = exp 



a(r)- 1 e e (r)d^ r - - / | C 7(r)- 1 ^ e (r)| 2 dr 



According to (HI), we easily get Eexp ^ J T |cr(r) 1 ^ e (r)| 2 dr < oo. By the Novikov condition 

and the Girsanov theorem, {Wt)o<t<T is a Brownian motion under the probability measure 
Q e := R e ¥. Then {X, X e ) is a coupling by change of measure with changed probability Q e . 
Since Rq = 1, by (39J Theorem 2.1], to obtain the desired integration by parts formula, it 
remains to confirm the following equality: in the sense of L X (P), 



— Re\e=0 - 

de 

Actually, noting that 

, ^ Re ~ 1 
limE— 

e^O e 



a(r)~ l \c H r^- H - rdb(r, -){X r 



T 



limE- 

limE- 

e->0 



Jo a(r)-^(r)d^ r - HI |q(r)-^(r)| 2 dr 
J T a(r)- 1 ee(r)d^ r ' 



and, moreover 



E 



+ 



T 



a{r)~ l C H r^ H - rdb(r, -)(X r ) 



-AWr 



< K 2 



E 



b{r,X?)-b(T,X r )-db(r,-)(X r We 



dr 



so the dominated convergence theorem implies the assertion. 

The second result in (1) follows by the given upper bounds on |cr(t) _1 | and \db(t, -)| and using 
the above integration by parts formula and the Young inequality (see, for instance, [2[ Lemma 
2.4]) 



\Pr{V y f)\ - a [Pr(f log f) - (P T f) (log P T f)} 
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< a log E exp 



ft 

< — log E exp 



- r^s)- 1 (c H s^- H -sdb{s,-){X s )) ^dW s 

J 

T 



ft- 



J \aisY 1 (c H s?- H -sdb(s,-)(X s 



)- 


2 

ds 


/ T 





■Prf 
Prf- 



Finally, (2) and (3) can be easily derived by applying [391 Proposition 2.3] and the second 
inequality in (1). The proof is complete. 

The shift Harnack inequality allows us to deduce the regularity for the law of the solution of 
(|3.ip . That is, we have the following result. 



Corollary 3.2 Suppose that the assumption (HI ) holds. Then, for any t > 0, the law of Xt is 
absolutely continuous with respect to the Lebesgue measure. 

Proof. Without lost of generality, we only consider the case t = T. Let 

pK 2 / n2 An ..is. , is2 



Ci 



r 2 



p - 1 V (2 - 2H)T 2H 5-2H ' 3 

The shift Harnack inequality stated in Theorem 13.11 implies, for any non-negative / S SS\, 

{P T f(x)fe~ c ^\ 2 < (P T {f(y + (x). 

Let A be a Lebesgue- null set, by applying the above inequality to / = L4 and noting the 
invariance property under shift for the Lebesgue measure, we have 



(P T Ia{x)Y [ e~ c ^ 2 dy<0, 

JR 



which implies the desired result. 

4 Bismut derivative formula 

In this section, we shall adopt the techniques of the Malliavin calculus to investigate the Bismut 
derivative formula and Harnack type inequality for Pt associated with (|3.ip . To this end, we 
make the following assumption: (H2) 

(i) there exist x G R and p > 2 such that b(-,x ) G L p ([0,T];R); 

(ii) b is differentiable w.r.t. the space variable such that db(t, ■) is uniformly continuous uni- 
formly w.r.t. the time variable t and moreover, 

\db(t,-){x)\ < K 3 ,K 5 < la(t)" 1 ! < K 4 , Mt G [0,T],x G R, 

where K^,K 4 and are positive constants. 

Main result reads as follows. 

Theorem 4.1 Assume that (H2) holds. Then, for all x,y G R and f G C£(R), 



V y P T f(x) = E 



f{X%)j\{s)- 1 (U + j\ H {s,r)u'{r)dr\ db(s, -)(X*) - u'(s)\ ydW s 



where u G C X (%T\;'&) such that 1 + K H (T,r)u'(r)dr = 0,X X is the solution of d3T 
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The proof of this theorem is based in the following lemmas and proposition. 

We first recall a result from [61 Theorem 4.1 and Theorem 4.2], in which the existence of 
Malliavin directional derivative is discussed. 

Lemma 4.2 Let b and a be continuously differentiable w.r.t. their space variable, with bounded 
derivative; assume further that, there exist xo G R and p > 2 such that b(-,xo) G L p ([0,T];W) 
and a is bounded. Then, for any £ G H, {{DnXf , C)n)t&[o,T] exists and is the unique solution to 
the equation 

Y t = (K H (K H (t,-)a(-,X x )),On + I K H (t,s)db(s,-)(Xf)Y s ds+ f K H (t, s)8a(s, -)(X*)Y s dW s , 

Jo Jo 

where X x is the solution of ()2.3I) . 

Following the same method presented in [61 Theorem 3.3], we can show that the solution of 
(|2.3D depends continuously on the initial condition in the sense specified below. 

Lemma 4.3 Assume b and a are Lipschitz continuous for the second variable uniformly w.r.t. 
their first variable, and there exist xo,yo G E such that b(-,xo) £ L+([0, T]; K), cr(-, yo) £ 
L^ 2 ([0,T];M). Denote by X x and X v the solution of ([23]) with initial condition x and y 
respectively. Then, for any p G Ah, there exists constant L p > such that 

sup E\Xf - Xf\ p < L p \x - y\ p . 
te[o,T] 

Remark 4.4 If we consider the case p = 2, Lemma \4.3\ reduces to {6[ Theorem 3.3]. 

Next we will concern the existence of the derivative process w.r.t. the initial data. 

Proposition 4.5 Suppose that b and a are both differentiable w.r.t. their second variables 
such that db(t, •) and da(t, •) are bounded and uniformly continuous uniformly w.r.t. their first 
variable t. Then, for each y G M, (VyX x )o<t<T exists and is the unique solution to the equation 

Y t = y+ f K H (t, s)db(s, -)(X*)Y s ds + f K H (t, s)da{s, -){X x s )Y s o\W s , 
Jo Jo 

where X x is the solution of ()2.3|) . 

Proof. Using the Picard iteration argument introduced in (6l Theorem 3.1], we can easily show 
that the above equation has a unique solution {Yt)te[o,T} an d moreover, sup ig [ 0T i E|Y^| P < oo 
holds for any p G Ah- 

For e > 0, let Z\ = X^ +ey — Xf — eY t , t G [0, T]. To complete the proof, it suffices to prove 

\Z e \ 2 

limEMj 1 - = 0, Vi G \0,T}. 
To this end, we see that, for any t G [0, T], 

Z\ = f K H (t,s){b(s,X x+ ^)-b(s,Xf)-edb{s,-)(Xf)Y s )ds 
Jo 
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+ f K H (t,s) (a(s,X*+ e y)-a(s,X*)-eda(s,-)(X*)Y s )dW s . 
Jo 

Therefore, by the Holder inequality and the Burkholder-Davis-Gundy inequality, there is some 
constant C 2 such that 

E|Z t f < 2T F K 2 H (t,s)\b(s,Xf +e v)-b(s,X*)-edb(s,-)(Xf)Y s \ 2 d S 
t 

+ 2C 2 [ K H (t,s)\a{s,X* +ey )-a(s,X*)-eda(s,-)(X*)Y s \ 2 ds 
Jo 

=: 2T [ K 2 H {t,s)J l {s) 2 ds + 2C 2 [ K 2 H (t, s)J 2 (s) 2 ds. 
Jo Jo 

Next we are to estimate Ji(s) and J 2 (s). Let us define, for each 5 > 0, 



(4.1) 



a(6) = sup sup {\db(s,-)(x) - db(s, -)(y)| + \da(s,-)(x) - da(s,-)(y)\) . 

\x-y\<8 se[0,T] 

It is clear from the assumptions on the coefficients b and a that a(oo) < 00 and a(5) 1 as 
5 4 0. As a consequence, we derive that, 

5 2 a 2 (5) = 5 2 a 2 (5)l {5<v - e} + S 2 a 2 (8)l {5>v - e} < 6 2 a 2 (V~e) + ^S^, 

e 2 

where g is chosen such that 2 < q < jj-. 
Note that, by the mean value theorem, we get 

Us) = I (db(s, -)(Ci) - 0&(s, -)PC)) (*? +ew - X?) + 56(«, -)(X*)Z% 

and 

J 2 ( S ) = I (da(s, -)(C 2 ) - -)(X*)) (X^y - X*) + &t( S , -){X*)Z% 

where Ci = 0*Xf + (1 - 6i)Xf +ty , Q % 6 (0, 1), i = 1, 2. 
Hence, we conclude that 



Ji(s) 2 + J 2 (s) 2 < (J^s) + J 2 (s)) 2 < 2a 2 (\X^ - Xf\)\X* +ey - X x s \ 2 + 2M\Z\\ 2 

< 2a 2 (V~e)\X^y - X*\ 2 + ^1\X^ - X^ + 2M|Z|| 2 , (4.2) 



e 2 

\2 



where M = sup se[0jT] (|<%(s, -)| + \da(s, -)\Y 
Now we turn to the estimate of E|Z|| 2 . Substituting (|4.2p into (|4.ip and noting q £ Ah, we 
have, by Lemma 14.31 



E|Z|| 2 < 2{TVC 2 )T 2H (2L 2 \y\ 2 a 2 {^7)e 2 + L q \y\ q a 2 [oo)e2 +1 ^ + A(T V C 2 )M J* K 2 H (t, s)E\ZI\ 2 ds 

= :C(e) + C 3 I K H (t,s)E\Z e s \ 2 ds. 
Jo 

We set K 2 (t,s) = K 2 H (t, s), K 2 +1 (t, s) = f* K 2 (t,r)K 2 (r, s)dr,Vs,t G [0,T],n > 1, and identify 
the operator K 2 with its kernel, K 2 := 1. 
Then by induction, we deduce that 

E|Z t f < C(e)(l + C 3 (K 2 l)(t)) + C 2 I f K 2 H (t,s)K 2 H (s,r)E\Ztfdrds 

Jo Jo 
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< ■■■<C(e)^Cl(Kfl)(t) + C^ 1 sup nZ^(Kl +1 l)(t). 
i=o ue [°> T ] 

Recall that [6j Lemma 3.3] states that X]£o su Po<t<T(^l)(0- 2;i < 00 ) ^ z e Therefore, 
letting n — >• 00, it follows that 

00 

E|Z t f < C(e)£C5 sup 
<=o * e [°' T l 



C(6) 

Now we are in position to prove Theorem 14.11 



Observing that lim — = 0, the proof is finished. 

€->0 e 



Proof of Theorem \4-l\ Note that if there exists £ G Dom<5f/ such that 

(D H X%,0 H = V y X%, a.s., (4.3) 

then for each / G CjfR), 

V y P T f(x) = V,E/(Af) = EV v /(Af)=E(/(Xf)V,Xf) 
= E(/'(X|)(£) H Xf,0«) =E((D H f(XZ),O n ). 

Applying the integration by parts formula for Dh, i.e. the definition of 5h, we get 



V y P T f{x) = E(f(Xft5 H £) = E (f(X%) £&HB* 



Furthermore, if £ G L 2 a {[0,T] x Q), then V y P T f(x) = E ) £ s dW s 

Based on the analysis above, we know that, to complete the proof, it suffices to find a £ = KhS, 

such that £ G L 2 a ([0,T] x fi) and JOJ holds. 

Let 

6 = ^a)" 1 ( (l + ^ r>' (r)dr) .)(Xf) - u'(s)\ y, 

where u G C 1 ([0, T];R) such that 1 + J Q T KH(T,r)u' (r)dr = 0. Obviously, £ constructed above 
is in L^([0,T] x 0). Next consider the following equation 

Z t = y+ f ' K H (t, S )db(s,-)(X*)Z s ds- f K H (t,s)a(s)tds. (4.4) 
jo Jo 

By the assumption, it is clear that (|4.4p has a unique solution Z. On one hand, observe that 
Yt := (1 + Jq Kff(t,r)u' (r)dr)y solves (|4.4p . On the other hand, since da(s,-) = 0,Vs G [0,T], 
Lemma [4.21 together with Proposition 14.51 implies that (V y X'f — (DhX£, £)H)te[o,T] is also a 
solution of (g3D. As a consequence, Y t = V y X? - (D H X£,£) H ,\/t G [0,T] holds. Due to F T = 0, 
it follows that {DjjX^,^)^ = X7 y X^. Therefore, the proof is complete. 

Remark 4.6 If we take u'(t) = — Hfrtz~ H , then by the proof of Theorem 13.11 we know that 



1 + Jq Kh(T, r)u'(r)dr = holds and the result of Theorem \4.1\ can be expressed as 



VyP T f(x) = E 



/(Xf ) £ o-(s)- 1 ((T - s)db(s, + |dW s 



In particular, when H = ^, we obtain a version of relation above that is an extension of \38\ 
Theorem 3.1], in which the coupling argument is used. 
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Next we will state some applications of the derivative formula obtained above. More precisely, 
explicit gradient estimate, Harnack inequality and log-Harnack inequality are presented. That 
is 



r 2 



Corollary 4.7 Assume that (H2) holds and set C(T,K 3 ,K 4 ,H) = 2K% , 

(1) For any f € 3S h {R), we get 

\V y P T f(x)\ 2 < C(T,K 3 ,K 4 ,H)\y\ 2 P T f 2 (x), 
i.e., \V v Pt f (x)\ is bounded above by f. Moreover, for all 8 > and positive f G 



(2-2H)T 2H J - 



\V y P T f(x)\ < 5 [P T (f log f) - (P T f) (log P T f)) + 



C(T,K 3 ,K 4 ,H). |2 



y\ 2 P T f. (4.5) 



(2) For any non-negative f £ i^,(R) and p > 1, the following Harnack inequality holds: 

P 



(PTf(x)f < P T f p (y)ew 



-C(T,K 3 ,K 4 ,H)\x-y\ 



(4.6) 



(4.7) 



p — 1 

As a consequence, the log-Harnack inequality 

P T (logf)(x) < log P T f(y) + C(T,K 3 ,K 4 ,H)\x-y\ 2 
holds for any positive f £ &b(M), and Pt is strong Feller, i.e. for each i£M, 

lim P T f(y) = P T f(x). 

(3) Let \i be Pt sub-invariant, i.e., fi is a probability measure on R such that J R Pr/d/U < 
J R fd/j, for all f € ^&(R), / > 0. Then the entropy-cost inequality 

»(P T f log P T f) < C(T, K 3 , K 4 , H)W${fn, v) 2 , f > 0, u.(f) = 1, (4.8) 
holds for the adjoint operator P T of Pt in L 2 (/j,), where d(x, y) = \x — y\. 

Proof. Let u'(t) = -^f-t^' H and define M T = / T a(s)~ l ((T - s)db(s, -)(X%) + C H s^- H ^j %dW s 
By the hypotheses on the coefficients, we derive that 



(M) 



T 



o-{ s y l ((T-s)db( s ,-)(x*) + c HS 



i- H \ y 



ds<C(T,K 3 ,K 4 ,H)\y\ 2 , 



where C(T, K 3 ,K 4 , H) = 2K\ (^p + {2 J» )T2H 
Hence, it follows from the Holder inequality that 

\V y P T f(x)\ 2 < E(M) T P T f 2 (x) < C(T,K 3 ,K 4 ,H)\y\ 2 P T f 2 (x). 

Combining the derivative formula with the Young inequality yield that, for any positive / E 
@ b (R) and 5 > 0, 



|V„JW(z)| < S[P T (f log f) - (P T f)(logP T f)]+5logEexp 



Prf- 



(4.< 
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Observe that 



Eexp 



M T 



< Eexp 



2{M) 



T 



< exp 



C(T,K 3 ,K 4 , H) ] 
8 2 



\y\ 



(4.10) 



Substituting (|4.10p into (f^T9|) implies (|4~5|) . In the sprit of [101 Corollary 1.2], (|4T6|) follows from 
(|4.5p . Since R is a length space, then according to [371 Proposition 2.2], (|4.6p implies (|4.7p . The 
strong Feller property follows from (|4.6p . due to the same proof of [8j Proposition 4.1]. Finally, 
(|4.8|) can be proved as the proof of [28^ Corollary 1.2] or [141 Corollary 3.6]. 



Remark 4.8 Making use of the Harnack inequality, one can compare the values of a refer- 
ence function at different points, while in the shift Harnack inequality presented in Theorem 
13.11 instead of initial points, a reference function is shifted. Besides, the (resp. shift) Harnack 
inequality allows us to compare the measure Pt(x, •) with some invariant probability measure as- 
sociated with a certain semigroup (resp. the Lebesgue measure), where Pt(x, •) is the transition 
probability for Pt- One can see J39} for more applications of the shift Harnack inequality. 



5 Transportation inequalities 

In this section we will discuss the Talagrand type transportation cost inequalities for the law of 
the solution of (12 . 3|) w.r.t. the uniform distance doo and the L 2 -distance cfo on the path space 
C([0,T];R). To the end, we introduce the following assumption: (H3) 

(i) there exists constant Kq(> 0) such that 

\b(t,x) - b(t,y)\ + \a(t,x) - a(t,y)\ < K 6 \x - y\, Vi E [0,T], x,yeR; 

(ii) \\c\\oo '■= sup sup \a(t, x)\ < oo. 

0<t<Tx&R 

Let us start by prove the following proposition which is crucial for the proof of Theorem 15.21 
below. 

Proposition 5.1 Let H > ^ and r be an {^t)stopping time. Assume that <j) is an adapted 
stochastic process satisfying E \<f>t\ p dt < oo for some p > 2. Then, there holds the maximal 
inequality 



E sup 

\0<t<TAr 



ft F\ rl At 

K H (t,s)<p(s)dW s j <C(p)E \cf> t \Pdt, 

where C(p) is a positive constant depending on p. 
Proof. Recall that, for H G (0, 1), Kn(t, s) is the kernel 

K H (t, s) = a H (t - s) H ~^ + a H Q - h) jT (r - sf-l (l - (^) dr, 



i _ ( 2HT(3/2-H) 

wnere a H — yr{H+i/2)r(2-2H) 
From this relation, we get 



1/2 



^-H)(rv»)»- 
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When H > \, the kernel K H (t, s) can reformulate as (for instance, see pQ and references therein) 



K H (t,s) = a H (^H - V\ S \~ H 1^ r H ~^(r - s) H ~*dr =: a H s^~ H 



IT -L , , IT d 

r 2 (r — s) 2 dr. 



Therefore, we have 



I K H (t,s)0(s)dW s = a H [ s^~ H [ r H ~^{r - s) H -^dr^(s)dW s . 

JO JO Js 

To exchange the integration of the right-hand side of the above expression, one need to show 
that the integrand fulfills the conditions of the stochastic Fubini theorem (see Theorem 4.18]). 
Actually, choosing e £ (0, i) such that H > and using the Holder inequality and the Young 
inequality, we obtain 



E s 1 - 2H (T-sf H - 3 ^(s)ds) dr 



< 



< 



< 



t \ 2 

-l+2e j 

r dr 



Jo Jo 



r — s 



2H 



T 2e \ 2 
"27J 

\4e(2H -e-1 



4H-2e-3 



drE / r^ Ze (f> z (r)dr 

T 







- 3 (P 2 (s)dsdr^j 



E / (j) z {r)dr 



o 



which is finite due to hypothesis on <f>. 
So, the stochastic Fubini theorem implies 



I K H (t,s)ct)(s)dW s = a H f ' r H ~h f t 
Jo Jo Jo 



{s)dW s dr. 



Taking # G (0, |) such that H > and applying the Holder inequality and the Young inequal- 
ity, we obtain 



( 



E sup 

\0<t<T/\T 
T 



r-t i-r P\ 

J r H ~^ J S ^- H (r - s) H ~l(f)(s)dW s dr J 



p-l r T 



E 



S h- H ( r - s) H -^(s)dW s 



l [o,TAr]{r) dr 



< 



< 



'E 



rAr 



S 2~ H (r - s) H ~2cp{s)dW s 



o 



dr 



p — 1 
Op 

p — 1 
0p 



p-l 



y>0p 



r(lr+ l_i_0 )p 



E 



^^-^(r-s^-^^dsVdr 



o 



p-l 



r — s 



fs) I/ e < r \I 



ds ) dr 



< 



< 



p 

1 \ P-l / 1 \2 ,2 , fTAT 

P~ l \ I \ T {H+\-e)\+{H-\) P -i^ I ,^ , OJ _ 



y (ff + i - 0)p + 2i? - 3 

which yields the desired result. 

We now prove the following main result in this section. 



•Rr, 
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Theorem 5.2 Let H > \. Assume (H3) and let F x be the law of the solution of (|2.3p with 
the initial point x on the path space C([0, T];R). Then, ¥ x satisfies the transportation cost 
inequalities on the metric space C([0,T];R). More precisely, 

(1) F x G T 2 (a(T,H)\d 00 ), where a(T, H) = 3(\\a\\ 00 T H ) 2 e 3K2 e T ( T2H + c ^ , 

(2) P x eT 2 (f3(T,H)\d 2 ), where (3(T,H) = 3(\\a\\^T H f^^^fl. 



Proof. Let Q be a probability measure on C([0,T];R) such that Q<P X . Clearly, to prove the 
desired result, we only need to consider the case H(Q|P a; ) < oo. The proof will divide into two 
steps. 

Step 1. The part follows the arguments of [II]. Let Q = J£(X)P. Note that 
(X)dP = / ^L( 7 )dP,( 7 ) = Q(C([0,T]; 



n air x ^c([o,T];TO^ dP, 



and 



that is, Q is a probability measure on (0, F) and H(Q|P X ) = M(Q|P). 

According to the proof of [HJ Theorem 5.6], there is a predictable process {ut)o<t<T such that 



^E f \u t \ 2 dt 



and the process 



o 



is a Brownian motion under Q, where Eq is the expectation taken for the probability measure 
Q. As a consequence, the process (Bf T )o<t<T defined by 



Bf = f K H (t, s)dW s = Bf - (K H u)(t) 
Jo 



is a Q- fractional Brownian motion associated with W. 
Step 2. From step 1, we can reformulate (|2.3p as 



X t = x+ [ K H (t,s)(b{s,X s ) + a(s,X s )u s )ds+ [ K H (t,s)a{s,X s )dW s . (5.1) 

JO JO 

Noting that, for any bounded measurable function F on C([0, T];R), 



Eq(F(X))=e(-^(X)F(X)) = / 



dl 



-( 1 )F( 7 )dP x ( 7 )=Q(F), 



it follows that the law of X under Q is Q. On the other hand, we consider the following equation 

Y t = x+ f K H (t, s)b(s, Y s )ds + f K H (t, s)a(s, Y s )dW s . (5.2) 
Jo Jo 
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As W is the Brownian motion under Q, we easily know that the law of Y. under Q is P^. 
Therefore, the law of (X, Y) under Q is a coupling of (Q, P x ) and moreover, we get 

W^°°(Q,P X ) 2 < E d oo (X,y) 2 = E ( sup \X t - Y t \ 2 ) , 

\0<*<T J 

w^(q,f x ) 2 < Eqc1 2 (x,y) 2 = e \x t - Y t \ 2 at^j . 

Combining (|5.ip with (|5.2h . we have 

X t -Y t = f K H {t, S )(b{s,X s )-b{s,Y s ))ds+ f K H {t,s)a(s,X s )u s ds 
Jo Jo 

+ f K H (t,s)(a( S ,X s )-a(s,Y s ))dW s . 
Jo 

Now, for n E N, define the stopping time 

r n := mf{t > 0, \X t - Y t \ > n}. 

Obviously, r n f oo as n goes to oo. Applying Proposition 15.11 and the Holder inequality, we 
derive that 

Eq(su P \(X t -Y t )l {t<Tn} \ 2 ) <3Eq(su P f ATn \K H (t, S )(b( S ,X s )-b( S ,Y s ))\ds) 

\0<t<T J \0<t<TJo J 

(rtATn \ 2 

sup / \K H {t,s)a(s,X s )u s \ ds) 
o<t<TJo l 



+ 3Eq sup 

\ 0<t<T 



< 3K 2 Eq sup 

\0<t<T. 



kt<r n y [ K H (t,s)(a(s,X s )-a(s,Y s ))dW s 
Jo 

rtATn \ 2 

J \K H {t,s)\-\X s -Y s \ds\ 
/ rtAT n \ 2 

+ 3||o-||;^Eq I sup / \K H (t,s)u s \ds) 

\0<t<TJo J 
/ rTAr n 

+ 3C(2)E (/ \a(s,X s )-a(s,Y s )\ 2 ds 

° r T 

< SQlaWooT 11 ) 2 ^ / u 2 s ds 
Jo 

+ 3K|(T 2 ^ + C(2))Eq f y |X S - y s | 2 ds 



= SdlalUr^)^ / u 2 ds 

+ 3K 2 (T 2H + C(2))Eq^ \(X s -Y s )l {s < Tn} \ 2 ds 



u 2 ds 










sup 1 




\0<t<s 
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By the Gronwall inequality, we obtain 



Eq( sup \(X t -Y t )l {t < Tn} \ 2 ) < 3(\\a\\ 00 T H ) 2 exp[3KlT(I aH + C(2))]K ii f 

\0<t<T J Jo 



u s ds. 



The Fatou lemma leads to 



sup \X t - Y t \ 2 ) < Sdlalloor^) 2 exp[3^|r(T 2H + C(2))]Ea / u 2 s ds. 

<t<T I JO 



\0<t<T 

Hence, we deduce that 

W$°°(Q,F X ) 2 < 2a(T,H)W(Q\F x ) 

with a(T, H) = 3(|| a\\ ^T 11 ) 2 exp[3K 2 T(T 2H + C(2))]. 

For the metric g?2, using the above procedure, we also can prove 



E (|(X t -y t )I{ t <r ?l }| 2 ) < 3(\W\\ooT H ) 2 ^ E q u 2 ds 

+ 3K 2 (T 2H + C(2)) ^ (|(X s - F s )I {s < Tn} | 2 ) ds. 

The Gronwall inequality, together with the Fatou lemma, yields 

X t -Y t \ 2 ) ^(IMUT") 2 rexp[3^ 2 (r 2 ^ + C(2))(t- S )]EQU 2 d S . 

Jo 



E 

Thus it follows that 

W^Q.P*) 2 < E ( r \Xt ~ Y t \ 2 dt 



< 3(\\ -\\ 00 T h ) 2 Eq u 2 ^ exp[3K 2 (T 2H + C(2))(t - s)]dt) ds 

< 2f3(T,H)U(Q\F x )°, 

where P(T,H) = 3{\\a\\ 00 T H ) 2 ^^ 

{T 2H +c{2)) ' P ro °f lc complete. 
Remark 5.3 In general, f Q tATn Kjj(t A r n , s)f(X s , Y s )dW s does not make sense, which forces 



us to consider Eq (J(Xt — Y t )I{ t<Tn y\ J rather than Eq (|Xf ATn — Y t ^ Tn \ 2 ). Further reading on 
stochastic Volterra equation, one can see \44] and references therein. 
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